Intermittency in spherical Couette dynamos 
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We investigate dynamo action in three-dimensional numerical simulations of turbulent spherical 
Couette flows. Close to the onset of dynamo action, the magnetic field exhibits an intermittent 
behavior, characterized by a series of short bursts of the magnetic energy separated by low-energy 
phases. We show that this behavior corresponds to the so-called on-ofF intermittency. This behavior 
is here reported for dynamo action with realistic boundary conditions. We investigate the role of 
magnetic boundary conditions in this phenomenon. 



I. INTRODUCTION 

First suggested by Joseph Larmor in 1919, 
dynamo action, i.e. the self-amplification of 
a magnetic field by the flow of an electrically 
conducting fluid, is considered to be the main 
mechanism for the generation of magnetic fields 
in the universe for a variety of systems, includ- 
ing planets, stars, and galaxies [l|. Dynamo ac- 
tion is an instability by which a conducting fluid 
transfers part of its kinetic energy to magnetic 
energy. 

In experiments, it is rather difficult to achieve 
a regime of self-excited dynamo action. The 
low value of the magnetic Prandtl number of 
liquid metals requires the injection of a suffi- 
ciently high mechanical power, and thus gen- 
erates turbulent flows, before reaching the dy- 
namo threshold. Dynamo action was first ob- 
served experimentally only in 2001, in Karl- 
sruhe [3 and Riga Q , and then in 2007 with a 
von Karman swirling flow of liquid sodium [J| . 

In parallel with these approaches, numeri- 
cal simulations have been carried out to model 
either laboratory experiments or astrophysical 
systems, for which the spherical geometry is rel- 
evant. We investigate spherical Couette flow 
and focus on the characteristics of the magnetic 
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field close to the dynamo onset. We observe a 
series of short bursts of the magnetic energy sep- 
arated by low-energy phases. This intermittent 
behavior, also known as on-off intermittency or 
blowout bifurcation, is usually interpreted as 
the effect of a multiplicative noise acting on a 
bifurcating system [a, |6|. 

On-off intermittency has so far never been 
observed in dynamo experiments, except in 
the case of an externally amplified magnetic 
field [7|. In contrast, it has been reported in 
a small number of numerical simulations |8l-ll0|. 
all relying on a flow in a periodic geometry pro- 
duced by a periodic analytic forcing. Here we 
investigate the influence of a realistic choice of 
boundary conditions on this phenomenon. 



II. GOVERNING EQUATIONS 

The spherical Couette flow geometry consists 
of two concentric spheres in differential rotation: 
the outer sphere, of radius Tq , is rotating around 
the vertical axis e^ with an angular velocity fl, 
and the solid inner sphere, of radius rj, is ro- 
tating at velocity Q + AQ around an axis that 
can make an angle 9 with e^. The aspect ratio 
X = Ti/ro is set to 0.35 to mimic that of Earth's 
liquid core. The spherical shell in between the 
two spheres is filled with an incompressible con- 
ducting fluid of kinematic viscosity v, electrical 
conductivity a, and density p. Its magnetic per- 



meability /iq is that of vacuum. The magnetic 
diffusivity 77 is defined as rj = l/(/iocr). 

We describe the problem in the reference 
frame rotating with the outer sphere. This in- 
troduces two extra terms in the governing equa- 
tions: the CorioHs force and the centrifugal ac- 
celeration. The latter can be rewritten in the 
form iV(rj^s^), where s denotes the distance 
to the axis of rotation. This term is a gra- 
dient and can be added to the pressure term 
which acts as a Lagrange multiplier to enforce 
the solenoidal condition on the velocity field. To 
establish the set of equations for this system, 
we rely on the same non-dimensional form as 
in [il|: the velocity u is scaled by rjAfi, the 
magnetic field B by y^p/ionro (fi -I- AH) Afi, 
and the length scale by Tq. The Navier-Stokes 
equation governing the fluid velocity u then 
takes the form 

du , ^, 2 , , 1 ^^ 



1 



1/1 Re 



Re Re \E x 



(VxB)xB, (1) 



and the induction equation for the magnetic 
field B, 

^.Vx,„xB) + ji-V=B^ (2) 

Both fields are solenoidal 



VB = 0. 



(3) 



The dimensionless parameters are the Ekman 
number E = v/ip^r^^, the Reynolds number 
Re = (rorjAri) /u, the magnetic Prandtl num- 
ber Pm = V /rj, and the magnetic Reynolds num- 
ber Rm = RePni. The potential 11 includes 
all gradient terms (the pressure term as well as 
the centrifugal effect introduced above). The 
Reynolds number varies with the rotation rate 
of the inner sphere, while the Ekman number 
is inversely proportional to the rotation rate of 
the outer sphere. When the latter is at rest, 
the Ekman number tends toward infinity and 
the Coriolis term in the Navier-Stokes equation 
vanishes. In our simulations, the Ekman num- 
ber is set to 10~^. This moderate value yields 
a moderate computing time. 



We impose no slip boundary conditions for 
the velocity field on both spheres. Magnetic 
boundary conditions are of three types. The 
first one can only be applied to the inner sphere, 
as it implies a meshing of the bounding solid 
domain. The inner sphere can be a conductor 
with the same electric and magnetic properties 
as the fluid. In that case the magnetic diffu- 
sion equation is discretized and solved in the 
solid conductor (we refer to this set of boundary 
conditions as "conducting"). The outer sphere 
as well as the inner sphere can be electrical 
insulators. In that case the magnetic field is 
continuous across the boundary and matches a 
potential field, decaying away from the bound- 
ary. The spherical harmonic expansion allows 
an explicit and local expression for these bound- 
ary conditions (we refer to these boundary con- 
ditions as "insulating"). In addition, the use 
of high-magnetic-permeability boundary condi- 
tions may enhance dynamo action |12| . There- 
fore, we also used boundary conditions which 
enforce the magnetic field to be normal to the 
boundary. This is equivalent to assuming that 
the medium on the other side of the boundary 
has an infinitely larger permeability (we refer to 
these boundary conditions as "ferromagnetic"). 
The different configurations investigated in this 
study are summarized in Tabled 

Inner Sphere Outer Sphere 
B.C.I Conducting Insulating 
B.C. 2 Insulating Insulating 

B.C. 3 Ferromagnetic Ferromagnetic 

TABLE. I. The three different configurations of 
magnetic boundary conditions used in this study. 



We integrated our system with parody jl3i] . 
a parallel code which has been benchmarked 
against other international codes. The vec- 
tor fields are transformed into scalars using 
the poloidal-toroidal decomposition. This ex- 
pansion on a solenoidal basis enforces the con- 
straints ([3]). The equations are then discretized 
in the radial direction with a finite-difference 
scheme on a stretched grid. On each concentric 



sphere, variables are expanded using a spher- 
ical harmonic basis (i.e., generalized Legendre 
polynomials in latitude and a Fourier basis in 
longitude). The coefhcients of the expansion 
are identified with their degree I and order m. 
The simulations were performed using from 150 
to 216 points in the radial direction, and the 
spherical harmonic decomposition is truncated 
at (/max, Wmax) = (70, 20). We obscrvc for both 
spectra a decrease of more than two orders of 
magnitude over the range of I and m. This 
provides an empirical validation of convergence. 
We checked on a few critical cases that the re- 
sults are not affected when the resolution is in- 
creased to /max = 100 . 

Let us define the non-dimensional kinetic and 
magnetic energy densities as 



Eh 



Ek 



1 1 



1 

v. 



1 



u dx. 



(4) 



K Re VE 



-) / B^dx, (5) 



in which the unit of energy density is p {vi Ai7) . 
In the above expressions, Vs refers to the vol- 
ume of the spherical shell. In addition, we also 
investigate the symmetry of the flow and the 
symmetry of the magnetic field with respect to 
the equatorial plane. To that end, we define 
the contributions to the energy densities cor- 
responding to the symmetric and antisymmet- 
ric components of the velocity (respectively E^s 
and EkA) and magnetic field (respectively E^s 
and Eh a)- The symmetric and antisymmetric 
contributions to the kinetic energy density re- 
spectively correspond to the flows 

us(a;, y,z) = ^ [u{x, y, z) + \i{x, y, -z)] , (6) 



ua{x, y,z) = \ [u(x, y, z) - u(x, y, -z)] . (7) 

In contrast, the symmetries are reversed for the 
magnetic field. This comes from the fact that 
the magnetic field is a pseudovector (i.e, the curl 
of a vector). Then, 

Bs(x, y,z) = \ [B[x, y, z) - B(x, y, ~z)] , (8) 



Ba{x, y,z)^\ \B{x, y, z) + B(x, y, -z)\ . (9) 

According to our definition, the dipolar compo- 
nent is symmetric. 



III. DIRECT NUMERICAL 
SIMULATIONS 

As shown by [11|, contra-rotation is more ef- 
ficient than co-rotation for dynamo action. In 
order to introduce more control over the sys- 
tem, we let the angle d between the axes of ro- 
tation of both spheres take any value in [0,7r]. 
Contrary to our expectations, we do not signifi- 
cantly lower the dynamo threshold with the in- 
clination of the rotation axis of the inner sphere. 
In fact, for Q = 7r/2, the fluid is mainly in 
co-rotation with the outer sphere, and dragged 
only by a thin layer on the inner sphere, which is 
not sufficient to trigger dynamo action. In our 
parameter regime, the best configuration seems 
to remain = tt, when the two spheres are in 
contra-rotation. We therefore keep this param- 
eter fixed in the rest of the study. 



A. Role of boundary conditions 

Let us first investigate the dynamo transition 
in this setup at fixed magnetic Prandtl num- 
ber Pm = 0.2, using the Reynolds number as 
the controlling parameter. With a conducting 
inner sphere and an insulating outer sphere, 
we find a critical magnetic Reynolds number 
Rmc € [291.0; 292.0], which is in good agree- 
ment with [ll|. Above the dynamo onset, the 
magnetic field displays an intermittent behav- 
ior characterized by series of short bursts of the 
magnetic energy separated by low energy phases 
(see Fig. [T]). When the distance to the thresh- 
old increases, bursts become more and more fre- 
quent and eventually intermittency disappears. 

Changing the boundary conditions generally 
leads to different thresholds for dynamo ac- 
tion. Using ferromagnetic boundary conditions, 
we find a critical magnetic Reynolds number 
Rmc e [298.6; 300.0]. With insulating bound- 
ary conditions, the threshold becomes large and 
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FIG. 1. Time evolution of the magnetic energy in linear (left) and log scale (right) for increasing Reynolds 
numbers at Pm = 0.2, using a conducting inner sphere (B.C.I in Table|l|. At lower Reynolds numbers (top), 
we see in linear scale the characteristic intermittent bursting. Intermittency gradually disappears at higher 
Reynolds numbers and the field reaches saturation. 



involves larger numerical resolutions. In order 
to maintain the hydrodynamic Reynolds num- 
ber at values which involve a moderate resolu- 
tion, we therefore had to increase the magnetic 
Prandtl number from 0.2 to 0.4 . We then ob- 
tain the dynamo onset for Rmc € [530.0; 534.8] . 
We emphasize that we observe the same inter- 
mittent regime with all the above choices of 
boundary conditions as long as the magnetic 
Reynolds number is close enough to the onset 
of the instability. 

For all boundary conditions, we observe 
that the dominant mode is predominantly of 
quadrupolar symmetry [the larger poloidal and 
toroidal modes are the {I = 2,m = 0) and 
(l ^ l,m = 0) modes, respectively]. For these 
Reynolds numbers, the flow is predominantly 
equatorially symmetric {EkA ^ Eks)- 




B. Increasing the magnetic Prandtl 
number 



Having assessed that the intermittent behav- 
ior of the magnetic field near onset could be ob- 
served with three different sets of boundary con- 
ditions, we restrict here our attention to simula- 
tions with ferromagnetic boundary conditions. 

Figure [2] presents the results we obtain at 
Pm = 2 . Close to the threshold, the mag- 
netic field still exhibits intermittency, but the 
nature of the process has significantly changed. 
There is now a clear distinction between two 
different regimes: phases of dynamo activity 
separated by phases of pure exponential decay. 
Both seem to alternate randomly. When the dy- 
namo is active, the magnetic field still displays a 
quadrupolar symmetry. In contrast, we observe 
the emergence of an axial dipole during decay- 
ing phases. The change of the global symmetry 
of the field coincides with the change of slope 
in the decaying phases [see Fig. [2] (bottom) and 
Fig. [3]. This change of slope is associated with 
a slower decay of the dipolar component over 
the quadrupolar mode. 



FIG. 2. Time evolution of the magnetic energy in 
linear and log scales for Pm = 2 and Re — 1200, 
using ferromagnetic boundary conditions. Instead 
of bursts, we now observe phases of dynamo activity 
which seem to randomly alternate with phases of 
exponential decay. The latter are no longer chaotic 
and are instead characterized by two different decay 
rates. 



IV. DISCUSSION 

A. Canonical model for on-ofF 
intermittency 

The simplest model that exhibits on-ofF inter- 
mittency is [iJI 



X = [a + C{t)] X-X^ 



(10) 



where a is the distance to the threshold, and ( 
a Gaussian white noise of zero mean value and 
amplitude D, defined as {C{t)C{t')) =^ DS{t-t') 
where () indicates the average over realizations 
(ensemble average). In the absence of noise, the 
system undergoes a supercritical pitchfork bi- 
furcation at a = . If a is sufficiently small, the 
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FIG. 3. (Color online) Time evolution of the sym- 
metric (dashed line) and antisymmetric (solid line) 
part of the magnetic energy. We focus on one of the 
decaying phases presented in Fig. (2] The decaying 
phase is characterized by a change of the dominant 
symmetry, as we can see on the visualizations of 
the magnetic field lines. The color insets respec- 
tively correspond to snapshots in the quadrupolar 
phase (left) and dipolar phase (right). 



B. Predictions and results 

Figure [4] shows the PDFs of the magnetic 
energy for a set of simulations at different 
Reynolds numbers. At low Rm, the PDF is 













10 







Rin = . -500.0 








--- 


Rm = 304.0 


8 








Rm = 30r).0 


■J-, 

o 




'' ^\ 




Rm = GOO.O 


bO 6 

O 




1 ~^ 




4 






^^"^ 


7. 


I 




?Ni ■ 




-11- 


10 -9 -8 


-7 


-6 -3 -4 -3 



logio Eb 

FIG. 4. Estimates of the probability density func- 
tions. Statistics are done from time series of 
the magnetic energy obtained with ferromagnetic 
boundary conditions, for Pm = 0.2. 



fluctuations lead to on-off intermittency, with 
bursts {a+( > 0) followed by decays (a+C < 0). 
During the off phases, one can neglect nonlin- 
earities and write Y = a+({t), with Y ~ \n{X). 
Thus, ln(X) should follow a random walk, with 
a small positive bias. Since solutions of Eq. ([T0| 
mimic solutions of the magnetohydrodynamics 
equations we observe in Fig.[Tl we further inves- 
tigate some properties of the model, (i) Equa- 
tion ([To]) leads to a stationary probability den- 
sity function (PDF) of the form |15| 



P(X) (X X(2-/^)-ie-^'/^ 



(11) 



which diverges at the origin for ^ .s = 
{2a/ D) — 1 < 1 . (ii) In addition, all the mo- 
ments of X must follow a linear scaling with 
a. (iii) Finally, another characteristic of this 
model is that the distribution of the duration of 
the off phases Toff follows a power law behav- 
ior, P{Toe) ~ T~g°', with a ^ -3/2. To com- 
pare these predictions to our results, we rely 
as in |10| on the magnetic energy density as a 
global measure of the magnetic field strength. 



characterized by a linear scaling on a log-log 
plot. The cutoff at low energies is not pre- 
dicted by the theory, which considers the limit 
i?b — >■ . For Rm > 310 , the magnetic energy 
fluctuates around a mean value and the PDF no 
longer scales as a power law. We see in Fig. [S] 
that the coefHcient s is proportional to the dis- 
tance to the threshold. Examples of the fit of 
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FIG. 5. Fit of the coefRcient s = {2a/ D) - 1, taking 
into account the linear domain of the PDFs in the 
intermittent regime only. Statistics are done from 
time series of the magnetic energy obtained with 
ferromagnetic boundary conditions, for Pm = 0.2. 
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the exponent s are presented in Fig. [Sj The 
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FIG. 6. (Color online) Examples of the fit (dashed 
lines) of the probability density functions (solid and 
dotted lines). Statistics are done from time series 
of the magnetic energy obtained with ferromagnetic 
boundary conditions, for Pm — 0.2. 

values of the coefficient are mainly affected by 
the range over which the data are fitted. Thus 
we select a range as large as possible. We then 
randomly sample this range with half-size sub- 
intervals. We then compute the mean slope 
and its standard deviation (represented in Fig.[S] 
with error bars). 

We then investigate the linearity of the mo- 
ments. Figure [7] shows our results for the first 
and second moments of the magnetic energy. 
We see that the mean magnetic energy grows 
linearly as a function of the magnetic Reynolds 
number. The second moments seem to follow 
the same linear trend, but only at the lower 
values of the magnetic Reynolds number. De- 
viations at larger values of Rm are expected, 
as this description is only valid in the limit 
Rm — >■ Rm^ . The duration tj of the time se- 
ries used to compute these values ranges from 
3.2 10^ to 1.410^ U.T. (values are presented 
in Table |ll| . These integration times are quite 
significant for a fully three dimensional set of 
partial differential equations but are necessarily 
short compared to the ones usually used with 
simplified models such as Eq. (fTO| . To quan- 
tify the uncertainty associated with the moment 
values, we sampled the integration time with 
sub-intervals randomly chosen. We then com- 
puted the moments on the full interval (symbols 



in Fig. [7]) and the standard deviation on the 
sub- intervals (reported as error bars) . The sub- 
intervals can be set from tj/A to tj/lQ without 
affecting these estimates. 



3.0 
2.5 
2.0 
1.5 
1.0 
0.5 
0.0 



300 302 



i' 



304 306 308 310 

Rm 



FIG. 7. (Color online) Moments of the magnetic 
energy as function of the magnetic Reynolds num- 
ber. Black squares represent the mean. The dashed 
line fits these data points with an error about 10% 
on the slope coefficient. The second moments (red 
circles) have been rescaled to match with the mean 
at Rm = 300 . 



Rm ti 


f//Rm 


300 1.4010* 


46.7 


302 9.0610^ 


30.0 


304 9.05 10=* 


29.8 


305 3.3310=* 


10.9 


306 7.58 10=* 


24.8 


308 3.20 10=* 


10.4 


310 5.4610=* 


17.6 



TABLE. II. Duration of the time series used to com- 
pute the moments of the magnetic energy. The in- 
tegration time ti is presented in units of (xAJl)" 
(ti) and in units of the magnetic diffusion time ra/rj 
{ti /Rm). 



Finally, we also tested the distribution of the 
duration time of the off phases. A definitive 
validation would require longer simulations, in 
order to have a significant number of off phases. 
For this reason, we can not rely on the simula- 
tions immediately above the threshold. Despite 



these short-comings, an iUustrative case is pre- 
sented in Fig. [8l Numerical values are given in 

Table nni 
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TABLE. III. Estimate of the exponent a of the PDF 
of the duration of the off phases for different thresh- 
old values and different ranges over which the fit is 
done. The standard error on the estimate of a is 
about 1%. Range values correspond to logj^o '^oS N 
axis in Fig. [S{b)[. 
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FIG. 8. Distribution of the duration of the off 
phases for a simulation performed at Pm ~ 0.2 and 
Re = 1470, with conducting boundary conditions. 
We define an off phase by a magnetic energy below 
a threshold value (indicated here by the horizontal 
black line) . The PDF (b) decay is close to a power 
law with the expected —3/2 value for the exponent 
(dashed line). 

To conclude, we emphasize that the predic- 
tions of the model are consistent with the three- 
dimensional simulations, and thus confirm the 
on-off hypothesis for the observed intermittency 
at low magnetic Prandtl number. 



Simulations at higher magnetic Prandtl 
number 



amining the dynamics of the flow. Indeed, we 
also carried out purely hydrodynamic simula- 
tions at Re = 1200 and observed intermittent 
transitions between two states. This kind of 
intermittent behavior of the flow was not re- 
ported in IllL but has been observed exper- 
imentally |16| . One state is characterized by 
larger fluctuations of the energy as we can see 
in Fig.[9l In addition, the analysis of the energy 




FIG. 9. Time evolution of the symmetric part of 
the kinetic energy Eks at Re = 1200 in a purely 
hydrodynamic simulation. Intermittent transitions 
between "laminar" and more turbulent phases are 
clearly visible. 



The simulations we performed at P,„ — 2 
exhibit a peculiar behavior of the magnetic 
field. This can be better understood by ex- 



spectra reveals that the ttt, = 3 modes domi- 
nate over the m — 2 modes during the "lami- 
nar" phases, whereas both are of the same order 



during the "turbulent" phases. Duration of the 
"turbulent" phases tends to increase gradually 
with the increase of the Reynolds number, so 
that the intermittent behavior of the flow even- 
tually disappears and is thus no longer present 
in the simulations at higher Reynolds number in 
which we have identified on-off intermittency. 
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FIG. 10. (Color online) Time evolution at Re = 
1200 and Pm = 2 of the total magnetic energy Ei, 
(dashed green line) and the kinetic energies E2 (dot- 
ted black line) and E3 (solid red line) for the m = 2 
and m = 3 modes, respectively. When the latter 
becomes larger than the former, the dynamo is no 
longer sustained and the magnetic energy exponen- 
tially decays. 



V. CONCLUSION 

Despite the fact that on-off intermittency has 
so far never been observed in dynamo exper- 
iments, we showed that the phenomenon can 
appear in numerical simulations of dynamo ac- 
tion using realistic boundary conditions. We 
identified in several cases the predicted behav- 
ior of the PDF of the magnetic energy, lin- 
ear scaling of the moments, and distribution 
of the duration of the off phases. In addi- 
tion, we tested these properties for three differ- 
ent boundary conditions (conducting inner core 
with insulating outer sphere, insulating or fer- 
romagnetic spheres). Finally, we pointed out 
a different kind of intermittency due to hy- 
drodynamic transitions that appears at lower 
Reynolds numbers. 

To explain the absence of on-off intermittency 
in the experiments, several reasons have already 
been invoked [131 • One explanation could be 
the imperfectness of the bifurcation (due for in- 
stance to Earth's ambient magnetic field). Since 
it has been shown that low-frequency noise con- 
trols on-off intermittency |14| . another possi- 
ble explanation could be that the low-frequency 
fluctuations are too small. However, the lack 
of experimental observations of on-off intermit- 
tency remains an open question and needs fur- 
ther investigations. 



Dynamo action is inhibited during the "lami- 
nar" phases (when the m = 3 modes dominate), 
which highlights the mechanism which leads to 
the peculiar behavior of the magnetic fleld, as 
we can see in Fig. [TOl In contrast, "turbulent" 
phases favor dynamo action, and one must wait 
a change in the flow to see the restart of dy- 
namo action after a phase of decay. Moreover, 
in a full magnetohydrodynamics simulation, we 
can artificially suppress the to = 3 modes of the 
velocity field by setting them equal to zero at 
each time step. We check that it is sufficient 
to suppress intermittency of the flow. Then, 
we observe that the phases of exponential de- 
cay are also suppressed and the dynamo is no 
longer intermittent. 
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